Sensitivity of Quantum Motion for Classically Chaotic Systems 
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importance which has been on the focus of many stud- 
ies p|-|ï^] is the so-called fidelity /(í), which measures 
the accuracy to which a quantum state can be recovered 
by inverting, at time t, the dynamics with a perturbed 
Hamiltonian: 



f(t) = Me* Mt e 
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We discuss the behavior of fidelity for a classically chaotic quantum system in the metallic regime. 
We show the existence of a critical value of the perturbation below which the exponential decay of 
fidelity is determined by the width of the Breit-Wigner distribution and above which the quantum 
decay follows the classical one which is ruled by the Lyapunov exponent. The independence of the 
decay rate from the perturbation strength derives from the similarity of the quantum and classical 
relaxation process inside the Heisenberg time scale. 

PACS numbers: 05.45-Mt, 05.45.Pq, 03.65.Sq 

Quantum chaos namely the attempt to understand 
classical dynamical chaos in terms of quantum mechanics 
has lead to a much better understanding of some prop- 
erties of quantum motion which go beyond simple in- 
tegrable models and perturbative treatments. A simple 
property of quantum conservative Hamiltonian systems 
with a finite number of partides, namely discrete spec- 
trum, has been at the origin of some difficultics. In- 
deed in the classical ergodic theory discrete spectrum to- 
gether with linear local instability of motion is a typical 
feature of integrable systems while chaotic systems are 
characterized by continuous spectrum and exponential 
local instability. This fact has questioned the possibility 
of dynamical chaos in quantum mechanics. On the other 
hand the correspondence principle requires transition to 
classical mechanics of all properties including dynamical 
chaos. As discussed in several occasions Q this apparent 
contradiction is resolved by taking into account that a 
sharp distinction between discrete and continuous spec- 
trum becomes meaningful only in the limit t — > oo. For 
finite times, there exist different time scales below which 
the quantum motion can display chaotic properties like 
the corresponding classical one. These time scales tend 
to infinite as the effective Planck constant ft e s — > 0. Two 
time scales are of particular importance: the random 
or the Ehrenfest time scale t r and the relaxation or the 
Heisenberg time scale tu- For t < t r the quantum motion 
is exponentially unstable like the classical one while the 
quantum relaxation process takes place during the time 
t < tfí. Since typically t r << tu, the quantum relax- 
ation process takes place in the absence of exponential 
instability. A clear illustration of this peculiar feature of 
quantum motion is shown in 0. It should be remarked 
that this lack of exponential instability does not prevent 
exponential decay of dynamical quantities like correlation 
functions or survival probability H . 

Recently the problem of the stability of quantum mo- 
tion has attracted a great interest, also in relation to 
the field of quantum computation. A quantity of central 
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Here ip is the initial state which evolves for a time t with 
the Hamiltonian Hq while H = Hq + V is the perturbed 
Hamiltonian. The analysis of this quantity has shown 
that, under some restrictions, the decay of f(t) is ex- 
ponential with a rate given by the classical Lyapunov 
exponent ||. This result appears to be consistent with 
recent experiments on the polarization echocs in nuclear 
magnètic resonance |Ï3| and with numerical computa- 
tions Q. More recent papers have contributed to clarify 
different complementary aspects of the problem |Mp[-p|, 
including the relation with the local density of states |Ï2] 
and the use of semiclassical approach (|] . The analysis of 
this quantity has some delicate aspects concerning some 
attempts to characterize quantum chaos via the classical 
Lyapunov exponent and the role of the above mentioncd 
time scales. It is therefore highly desirable to have very 
accurate numerical results and to this end it is necessary 
to consider simple systems which display the genèric fea- 
tures of classical and quantum chaotic systems and which 
can be easily treated numerically. 

In this paper we consider the behavior of fidelity for 
a classically chaotic system, in the delocalized regime of 
quantum ergodicity, in which the wave functions have 
a complex pattern which can be described within the 
framework of random matrix theory. We show that the 
type of decay and its rate depend on the strength of the 
perturbation. In particular, above a critical border, the 
quantum decay mimics the classical one and therefore, 
up to the relaxation time scale, it follows the exponential 
classical decay, which in our case is ruled by the Lya- 
punov exponent. The independence of the decay rate on 
the perturbation, which takes place in this regime, simply 
refleets the properties of the underlying classical motion. 
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We consider the classical sawtooth map: 

n = n+k (9-ir), 6 = 0+Tn, (2) 

where (n, 8) are conjugated action-angle variables (0 < 
9 < 2tt), and the bars denote the variables after one map 
iteration. Introducing the rescaled momcntum variable 
p = Tn, one can see that the classical dynamics depends 
only on the single parameter Kq = k^T . The map (|^) 
can be studied on the cylinder [p S (— oo,+oo)], which 
can also be closed to form a torus of length 2irL, where 
L is an integer. For Kq > the motion is completely 
chaotic and diffusive, with Lyapunov exponent given by 
A = ln[(2 + Kq + ((2 + Kq) 2 - 4) 1/2 )/2]. For Kq > 1, the 
diffusion coefficient is well approximated by the random 
phase approximation, D ps (tt 2 /3)Kq. 

The quantum evolution on one map iteration is de- 
scribed by a unitary operator Uq acting on the wave func- 
tion ip: 

^=Ü ^ = e -iTn 2 /2 e iko(Ò-*f/2^ ( 3 ) 

where n = —idjdB (we set % = 1). We take —N/2 < 
n < N/2, k = (K /2irL)N, T = 2tïL/N. The classical 
limit corresponds to N — > oo. We note that in this sim- 
ple quantum model one can observe important physical 
phenomena like dynamical localization and cantori local- 
ization H] . Our aim is to study the fidelity decay in the 
delocalized regime of quantum ergodicity. Moreover we 
will consider parameter vàlues for which there is no ini- 
tial transient diffusive behavior, which may considerably 
affect the decay of fidelity. 

In order to compute the fidelity we choose to per- 
turb our system by slightly varying the kicking strength, 
K = Kq + e, with e <C Kq. Correspondingly the per- 
turbed quantum kicking parameter is k = fco + a, with 
a = eN / (2nV). Since we want to compare classical and 
quantum evolution, we compute the classical "fidelity" 
f c {t) in the following way: we consider in the phase space 
a uniform density of points inside a strip of area A = 2nv 
(0 < e < 2?r, -v/2 <p< v/2). We then define f c (t) as 
the overlap of the initial area A with the area A' ob- 
tained by evolving A for t iterations of the map (2) and 
then reversing the evolution for t iterations with the per- 
turbed strength K = Kq + e. In practice, we follow the 
evolution of 10 6 trajectòries uniformly and randomly dis- 
tributed inside the area A and define the fidelity / c (í) as 
the percentage of orbits which return back to the area A 
at time t, after the above reversing procedure. The corre- 
sponding quantum initial condition is given by a uniform 
mixture of momentum states located inside the area A. 
We note that this choice, besides giving the correct clas- 
sical limit when N — > oo, introduces a convenient averag- 
ing procedure. Moreover, we have checked that the same 
fidelity decay rates are obtained if one starts from pure 
states, like momentum eigenstates or coherent states. 




FIG. 1. Decay of classical fidelity for the classical saw- 
tooth map with Kq = 1, L = 1, v = 27r/10 4 , and pertur- 
bation strength e = 10 -3 (circles), 10 -4 (squares), 10 -5 (di- 
amonds), 10 -6 (triangles), and 10 -7 (stars). The straight 
lines show the decay / c (í) oc exp(— Aí), with Lyapunov expo- 
nent A = 0.96. The dashed line indicates the saturation value 
/ c ,oo = v/(2ttL) — 10~ 4 . Here and in the following figures 
the logarithms are decimal. 
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FIG. 2. Decay of the fidelity for the quantum sawtooth 
map at Ko = 1, L = 1, e = 5 x 10~ 5 , N = 8192 (dot- 
ted line, a — 0.065), 16384 (dashed line, a = 0.13), and 
32768 (sòlid line, a = 0.26). The straight line gives the decay 
f(t) = exp(— Ft), with rate V = 2.2a 2 . As initial state we 
take a momentum eigenfunction with n = 0. 

The behavior of the classical fidelity is shown in Fig.l, 
for Kq — 1, L — 1, and different vàlues of the perturba- 
tion strength e. In this particular regime, characterized 
by (i) uniform local exponential instability and (ii) ab- 
sence of diffusive regime, the fidelity decay is ruled by the 
Lyapunov exponent A. The exponential decay starts after 
an initial transient time proportional to ln(v/e), which is 
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required to amplify the perturbation up to the scale v 
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FIG. 3. Classical and quantum fidelity decay for Ko = 1, 
L = 1, v = 27r/10 3 . Left curves: e = 1(T 3 , N = 16384 
(dashed line, a = 2.61), N = 131072 (sòlid line, a = 20.9), 
and classical decay (circles). Right curve: e = 10 -4 , 
N = 131072 (sòlid line, a = 2.09) and classical decay (tri- 
angles) . 
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FIG. 4. Rate 7 of the exponential decay for the quantum 
fidelity versus perturbation strength cr, for Ko = 1, iV = 2048 
(circles), 8192 (diamonds), and 65536 (squares), Ko = 2, 
N = 8192 (stars), K = 10, N = 8192 (triangles). The 
dashed line gives the decay rate F — 2.2a 2 . The sòlid lines 
show the Lyapunov decay, with rates A = 0.96 (at Ko = 1), 
1.32 (K = 2), and 2.48 (K = 10). 

The decay of the quantum hdelity is Gaussian below 
a perturbative border This border is given by the 

value of the perturbation at which the typical transition 
matrix element U between quasienergy eigenstates be- 
comes larger than the average levels spacing 1/p. For 



ergodic eigenfunctions, U ~ cr/VN, while the density 
of quasienergy states is given by p = N/2%. There- 
fore the perturbative border is given by a v w \j\fN. 
Above this border one typically expects an exponential 
decay of fidelity, with a rate T = 2npU 2 « cr 2 given by 
the width of the Breit-Wigner local density of states M . 
This theoretical prediction is confirmed in Fig.2, which 
shows the decay of quantum fidelity at e = 5 x 1CP 5 
and diffcrcnt N vàlues, with cr > a p . The nice scaling 
behavior of Fig.2 confirms the predicted exponential de- 
cay /(i) w exp(— Ca 2 i), with the numerically determined 
constant C w 2.2. 

On the other hand, as stated in the introduction, 
one expects that in the semiclassical regime the quan- 
tum motion mimics the classical one up to the relax- 
ation time scale which is determined by the density of 
quasienergy eigenstates which signihcantly contibute to 
the wave function dynamics. To this end it is neces- 
sary that the perturbation a is strong enough to allow 
the quantum motion to follow, on the average, the ini- 
tial classical decay. In our case this may happen if a is 
large enough to induce transitions at least between near- 
est neighbors momentum states, namely a > u c w 1. 
If a < a c , the quantum excitation is unable to follow 
the classical spreading of the initial state. One may also 
argue in a different way: since with our choice of param- 
eters we are in the metallic regime, all N quasienergy 
states are involved in the evolution of the unperturbed 
system. Then the cffcct of the perturbation on the quan- 
tum motion can imitate the corresponding classical one 
only if there are no quantum localization effects on the 
quasienergy states. This happens when the width of the 
local density of states becomes comparable to the band 
width, that is pT w N, which again gives the threshold 
value <7 C w 1. We remark that, as discussed in Jï^|, in 
the theory of Wigner band random matrices the Breit- 
Wigner regime corresponds to a sort of partial perturba- 
tive localization. The above theoretical estimate is well 
confirmed by our numcrical data presented in Figs.3,4. 
Fig. 3 shows that for cr > 1 the quantum fidelity follows 
closely the classical behavior, namely it decays exponen- 
tially with the classical rate given by the Lyapunov ex- 
ponent. Fig. 4 shows the decay rate 7 as a function of 
the perturbation strength a. It is clearly seen that for 
cr < 1 the decay rate is proportional to a 2 , that is to the 
width of the Breit-Wigner. Therefore a c w 1 is a critical 
value which separates two distinct regimes: a pure quan- 
tum perturbation dependent regime, and a semiclassical 
regime. We note that the perturbation a depends on 
both iV and e. For a > 1, the decay rate does not change 
by increasing N at fixed e, since by doing this we merely 
increase the Heisenberg time. On the other hand, if wc 
increase e at fixed N (provided that the perturbation re- 
mains classically small, i.e. e << Kq) the decay rate also 
does not change, since the exponential amplification of 
the perturbation is controlled by the parameter K w Kq. 
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In both cases the decay rate of fidelity is perturbation 
independent. This is a property of the classical motion 
which, in the semiclassical regime, is shared by quantum 
mechanics. However, we would like to stress that the 
decay of fidelity remains perturbation dependent, since 
the exponential decay starts after a time oc |lne| (see 
Figs.1,3). 
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FIG. 5. Fidelity decay in the diffusive regime with L = 50, 
Ko = 1, v = 7T, e = 0.05. The dashed curve gives the classical 
decay, the sòlid curves the quantum decay: from top to bot- 
tom N = 1024 («7 = 0.16), 2048 (cr = 0.33), 8192 (cr = 1.30), 
16384 (a = 2.61), and 32768 (cr = 5.22). The straight line 
indicates the decay f(t) oc l/V Dt. 

For the parameters vàlues of Figs.1-4, the decay of 
fidelity is exponentially fast and the saturation value 
foo = v/{2"ïïL) is reached on times much shorter than 
the Heisenberg time. In order to observe the effect of 
the Heisenberg time scale it is necessary to have a much 
slower decay of fidelity. In Fig.5 we take Ko = 1 and 
L = 50, so that we allow for a Gaussian diffusive pro- 
cess in momentum space. Because of this, during the 
diffusion time the fidelity decays in the classical case as 
\/\[~Dt Fig.5 shows that for a > a c as 1 the quan- 
tum decay follows the classical one for larger and larger 
times as N increases, in agreement with the correspon- 
dence principle. The asymptotic value is foo = vi/ (2ttL), 
where, according to the scaling theory of localization, 
/ = Ç/N = g(x), with x = k 2 /N |]|. Here £ is the ac- 
tual localization length of the "sample" of size N, while 
k 2 gives the localization length for the infinitc sample, up 
to a numerical constant of order 1. The scaling function 
g(x) is proportional to x for x -C 1 and saturates to 1 
for x ^> 1. The transition value x — 1 corresponds to 
N ps 10 5 . Moreover, the saturation value is approached 
after a relaxation time ti ~ £. We stress that in the case 
of Fig.|| the decay of fidelity is controlled by the diffusion 
coefficient and not by the Lyapunov exponent. The ob- 
servation of such regime represents a challenge for experi- 



ments like spin echoes. Further theoretical investigations 
are also desirable in order to understand more clearly the 
effect of classical diffusion and quantum localization on 
the behavior of fidelity. 

In summary, we have shown that the decaying behav- 
ior of fidelity in a classically chaotic system strongly de- 
pends on system parameters as well as on the perturba- 
tion strength. Nevertheless there is a regime in which 
the decay rate (exponential or power law) is perturba- 
tion independent: in this regime the quantum motion 
simply mimics the properties of the underlying classical 
dynamics. We emphasize that the quantum to classical 
correspondence of the average behavior is vàlid until the 
Heisenberg time scale, which is much longer than the 
Ehrenfcst time scale associated with the exponential in- 
stability of quantum motion. 
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